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ABSTRACT 

We present black hole mass measurements from kinematic modeling of high-spatial resolution 
integral field spectroscopy of the inner regions of 9 nearby (ultra-)luminous infrared galaxies in 
a variety of merger stages. These observations were taken with OSIRIS and laser guide star 
adaptive optics on the Keck I and Keck II telescopes, and reveal gas and stellar kinematics 
inside the spheres of influence of these supermassive black holes. We find that this sample of 
black holes are overmassive (~ 10 7 ~ 9 M 0 ) compared to the expected values based on black hole 
scaling relations, and suggest that the major epoch of black hole growth occurs in early stages 
of a merger, as opposed to during a final episode of quasar-mode feedback. The black hole 
masses presented are the dynamical masses enclosed in ~25pc, and could include gas which is 
gravitationally bound to the black hole but has not yet lost sufficient angular momentum to be 
accreted. If present, this gas could in principle eventually fuel AGN feedback or be itself blown 
out from the system. 

Subject headings: Galaxies: evolution - galaxies: kinematics and dynamics - galaxies: nuclei - galaxies: 
interactions 


1. Introduction 

Virtually every massive galaxy hosts a super- 
massive black hole in its core. These central black 
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holes exhibit tight correla tion s with p roperties of 
their host galaxies’ bul ges (IKormendv fc Gebhardt 
2001): bulge mass (IKormendv &; Richstonel 19951 


Marorrian et al .1 19981) . bulge luminosity (lMarami_^_Himt| 

2003), and bulge stellar velocity d i spersion (ITremaine et al. 


2002 ; 
20001 ) . 


Ferrarese fc Merritt] l2000t 1 Gebhardt et al 


Recently these correlations were re¬ 
explored using the wealth of new black hole mass 
me asurements (10 6 ~ 10 Mm) and gala xy para m eters 
by McConnell fc Mai ( 20131) and lGraham fc Scott 
(|2013l ). However, the discovery that t he Mbh — ct. 


relati on may evolve with redslrift (e.g. Zhang et al 


20121 ) has led to the suggestion that the most basic 
scalin g relation may instead be with total stella r 
mass ( Jahnke et al. 20091 : Cisternas et al. 2011a ): 
in this case, the evolution of the Mbh — v* rela¬ 
tion with redshift indicates the changing fractions 
of mass in galaxy bulges (which contribute to <r») 
versus disks (which do not). See IKormendv &; Ho 
( 20131) for a detailed review of our current under¬ 
standing of black hole scaling relations. 

The mechanism through which black hole 
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masses correlate with galaxy properties has been 
canonicall y associated wi t h gas -rich galaxy merg¬ 
ers (e.g. Hopkins et, al.1 l200(il) . Gravitational 
torques funnel the gas into their centers, trig¬ 
gering two phenomena: an intense burst of star 
formation to feed the bulge, and accretion of gas 
on to the black holes in the centers of each galaxy. 
It has been postulated that black hole growth 
can regulate this process through AGN feed back 
( Sprineel et al. 2005a : Hopkins fc Elvis \‘2() 1 Ol) via 
massive winds that evacuate the gas from the 
galaxy on short timescales, cutting off star forma¬ 
tion and future black hole growth. This sense of 
self- regula tion h as bee n con firmed observationally 
by iKauffmann fe Heckmanl ( 2009h , who find that 
the Eddington ratio of a sample of AGNs depends 
on the supply of cold gas in the galaxy. If there is 
plenty of cold gas, the accretion rate does not de¬ 
pend on the quantity of gas available; if the supply 
of cold gas is limited, the accretion rate depends 
on the rate at which stellar winds provide fuel 
for the AGN. Though the detailed mechanisms 
causing these correlations are still unconfirmed, 
star formation and black hole growth are fed by 
the same reservoir of inflowing gas; their growth 
histories are intertwined. It is likely that these 
two processes compete for fuel in a predictable 
fashion. 

To understand this interplay, it is critical to 
look at systems in the midst of this increased fu¬ 
eling. One set of such galaxies are gas-rich merg¬ 
ers, which tend to have extreme bursts of star 
formation and a higher incidence of AGN activity 


formatio n and a nigner incidence ot Au.N activity 
(e.g. Sanders et a ll 19881 S anders fc Mirabell 19961: 


Veilleux et al.1 2 )02t llshidal l2004 Ellison et al. 
20131 : Koss et al. 2013 ). During such a merger, 
does the black hole grow first, leaving the stars 
to slowly consume the remaining gas? Or is star 
formation quenched once the black hole reaches a 
bright quasar phase of extreme growth? The po¬ 
sition of a merger on black hole scaling relations 
would indicate the relative growth timescales, and 
confirm whether the putative quasar-mode feed¬ 
back occurs at the end of a merger (see Figure[l]for 
the schematic example for the Mbh—&* relation.). 
As star formation happens on few-kiloparsec spa¬ 
tial scales, accreting gas must lose less angular 
momentum to fuel a starburst than to feed a black 
hole. Subsequent models therefore suggest that a 
black hole would grow substantially only after star 


formation h as qu enched itself and the galaxy bulge 
is in place (Cep 120121 ), or at least that the peak 
black hole accretion time occurs l ater in a m erge r 
relative to peak star formation ( Hopkinsi 20121) . 
This scenario would predict that gas-rich mergers 
would fall below black hole scaling relations. 



a, (km/s) 


Fig. 1. — Mrh — cr* relation for isolated galaxies 
from McConnell fc M ;il ( 2013i ) (black) with three 
possible evolutionary tracks for merging galaxies 
overlaid. If the black hole grows first or more 
quickly than the galaxy bulge, mergers would lie 
above the relation (as shown by the gold arrows). 
If the black hole growth lags the bulge growth and 
is responsible for curtailing evolution (e.g. quench¬ 
ing through AGN feedback), mergers would lie be¬ 
low the relation (as shown by the purple arrows). 
If instead the black hole and the bulge grow in 
lockstep, the mergers would remain on the rela¬ 
tion (as shown by the blue arrow). 


Merger-driven galaxy evolution is not a com¬ 
plete explanation for all black hole growth, how¬ 
ever. Though some AGN st udies find a cor¬ 
relation with maj or mergers ( Koss et al. I 201(1 
Elli son et al ] 2011) , AGNs found with other selec¬ 
tion techniques and at different redshift ranges do 
not show a higher rate of merging than field galax¬ 
ies ( Cisternas et al.l 1201 lbl : iKocevski et al. 2012 ). 
Instead of mergers bringing in gas, some sys¬ 
tems probably are undergoing secular evolution, 
accreting their gas directly from the cold inter- 
galactic medium. Bar and spiral disk instabilities 
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are capable of dissipating sufficient angular mo- 
menturn from t his gas to fuel black hole growth 
( Hopkins fc Quataertll2010al 2011 1. The dividing 
line between these two processes is not yet well 
understood; here we focus on understanding the 
black hole growth due to major mergers. 

Nearby (Ultra-) Luminous InfraRed Galaxies 
((U)LIRGs: ISanders fe Mirabell 19961 ) are an ex¬ 
cellent sample with which to study such gas-rich 
mergers. These galaxies have infrared luminosi¬ 
ties upwards of 10 11 Lq (10 12 Tq for ULIRGs), 
generally caused by a starburst and/or an AGN 
heating up dust. Their infrared luminosities cor¬ 
relate with mer ger rate, star forma tion rate, and 


AGN fraction ( Ellison et a]. 20131 ). In fact, in 


the local universe, such strong infrared activity 
is alm ost ex clu sively triggered by major mergers: 
Veilleux et al. ( 20021) showed that in a complete 
sample of IRAS ULIRGs, 117 out of 118 galaxies 
are in the midst of strong tidal interactions. 

In a companion paper (Medline et al.l[2014 ), we 
have studied the kinematics of both gas and stars 
in the inner kiloparsecs of a sample of (U)LIRGs, 
finding that nuclear disks on scales of a few tens 
to hund reds of parsecs are common. In two other 
papers ( Medline et al. 2011 : U et al. IB , we 
demonstrated a technique that uses high spatial 
resolution integral field spectroscopy to measure 
black hole masses. With kinematic maps that re¬ 
solve inside the sphere of influence of a black hole, 
the complex and unrelaxed large-scale dynamics 
are less important, and black hole masses can be 
measured to within a factor of a few. We note 
that this technique measures the unresolved cen¬ 
tral mass, which includes both the black hole and 
its accretion disk, and in some cases may also in¬ 
clude a res ervoir of gas feeding the accretion disk. 
In U et al. ( 2013 ). the black hole mass in Mrk273 
N measured with this technique was consistent 
with the m easurement made b y OH maser kine¬ 
matics ([Klpckner & Baan 120041 ). 

Using a robust technique such as this to mea¬ 
sure black hole masses in these gas-rich mergers 
is important because more traditional methods of 
black hole mass measurements rely on assump¬ 
tions that aren’t valid in the case of galaxy merg¬ 
ers. Three-integral orbita l superposition mo dels 
(as in Giiltekin et al. 2009t Sioois et ~ahl 120091) are 


however, this approach requires a dynamically re¬ 
laxed system and is used therefore in isolated 
galaxies. Another successful black hole mass mea¬ 
surement technique is reverberation mapping (e.g. 
iDennev et al. 20091 and references therein), which 
measures the time lag between flux variations of 
the continuum and the lines in the broad line re¬ 
gion. Since the cores of (U)LIRGs are so dusty, 
the broad line regions are too obscured to view. 

Obtaining data at spatial resolutions suffi¬ 
ciently high to resolve inside the sphere of in¬ 
fluence of a supermassive black hole at the typi¬ 
cal redshifts of local (U)LIRGs (z < 0.1) requires 
adaptive optics systems, which are becoming avail¬ 
able at an increasing number of ground-based ob¬ 
servatories. Though the Hubble Space Telescope 
has excellent resolution in the visible bands, its 
relatively small mirror size limits the resolution 
at the longer wavelengths (> 2/rm) necessary to 
look through the dust in these galactic nuclei; 
even at the longest wavelength available to the 
Wide Field Camera 3 (H-band), these nuclei are 
still sometimes obscured. In order to achieve high 
spatial resolution in Tv-band, large ground-based 
telescopes have employed adaptive optics systems 
which measure turbulence in the Earth’s atmo¬ 
sphere and use a deformable mirror to correct 
for the resulting distortions. These distortion 
measurements require references, either a natu¬ 
ral guide star (NGS AO) or a laser guide star plus 
a fainter natural “tip-tilt” star (LGS AO). The 
addition of laser guide star adaptive optics has 
increased the area of the sky observable with this 
technique. 

Throughout this paper we have adopted a cos¬ 
mology of H 0 = 70kms~ 1 Mp c -1 , = 0.28, 


and Ua = 0.72 ( Hinshaw et al. 20091) . In (J2] we 


able to use large-field kinematics to separate differ¬ 
ent components to the mass profile of the galaxy; 


present our data and reduction techniques. In !}3] 
we briefly descri be the kinematic fittin g te chniques 
demon strated in lMedling et al.l (l201ll ) and II et al. 
d2013l) . In g] we present the black hole masses 
measured from several tracers and in fj5] compare 
them to black hole scaling relations. (J6] contains 
our conclusions. 

2. Observations 
2.1. The Sample 

We have selected 9 gas-rich merging galaxies 
in which to measure the black hole masses using 
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gas and stellar kinematics. These galaxies repre- 
sent the subset of m erging galaxies presented in 
Medline et al. ( 2014 1 for which high quality kine¬ 


matics data exist. The parent sample was drawn 
from the Great Observator ies All-Sky LIRG Sur¬ 
vey (GOALS: lArmus et al.ll2009h . which targeted 
about two hundred of the brightest infrared galax¬ 
ies in the sky (log(L/fl/L 0 ) > 11.0). Our targets 
were selected from that sample for the following 
criteria: available archival B- and /-band imaging 
from the Hubble Space Telescope Advanced Cam¬ 
era for Surveys; visible double or extended cores 
or clumps in the central five arcseconds (within 
our field of view); large-scale morphology consis¬ 
tent with that of a major merger; observable from 
Keck Observatory and with an appropriate guide 
star (< 18th magnitude in R, < 60" separation). 

Our sample comprises (U)LIRGs involved in 
major gas-rich galaxy mergers, with a bias towards 
later stages of merging. These galaxies are nearby, 
with redshifts < 0.05, which is required to achieve 
spatial resolutions of < 30 pc per pixel (< 50 pc 
per resolution element). 


2.2. The Data 


We obtained near-infrared integral field spec¬ 
troscopy of the central kiloparsec of 9 merging 
galaxies with OSIRIS, the OH-Suppressing In- 
fraRed Imaging Spectrograph ( Larkin et al.ll2006 ) 


on the W. M. Keck I (after August 2012) and 
II (before August 2012) 10-meter telescopes. All 
data presented here uses the 0.035 arcsec spaxel -1 
plate scale. This high spatial resolution is en- 


tern (IWizinowich et al. 

200(1 van Dam et al. 

2004; 

Wizinowich et al. 2006 

: van Dam et al. 20061). 


Our OSIRIS observations were typically com¬ 
prised of observing sets of object-sky-object; each 
exposure was ten minutes (5 minutes for Keck I ob¬ 
servations). Observations were taken either in the 
broad-band K filter (Kbb: 1.965-2.381 /im) or in a 
narrow band targeting specific lines (Hn4: 1.652- 
1.737 /im; Kn5: 2.292-2.408 /im). When possi¬ 
ble, the primary target lines were the CO (2-0) 
and (3-1) bandheads at 2.293 /im and 2.323 /im 
in K , which trace the kinematics of young stars. 
When available, we compare the stellar kinemat¬ 
ics to those of various emission lines: [Fe n] 1.644 
/im, Bry 2.16 /im, H 2 2.12 /inr. For some galaxies, 
stellar kinematics were not available because the 


CO bandheads are redshifted out of K , and we 
rely on the aforementioned emission lines for this 
analysis. Total exposure times, observed filters, 
and parameters of the observed galaxies are listed 


in Table fTl Some of these data were 1 

resented 

in previous papers ( 

Medling et al. 2011; 

lU et al. 

2013; Medling et al. 

2014), these are indicated in 


Tabled] as well. 


We reduced our Keck II observations using the 
OSIRIS Data Reduction Pipeline^] version 2.3 (us¬ 
ing the updated OSIRIS wavelength solution for 
data taken after October 2009) and our Keck I ob¬ 
servations using version 3. This pipeline includes 
modules to subtract sky frames, adjust channel 
levels, remove crosstalk, identify glitches, clean 
cosmic rays, extract a spectrum for each spatial 
pixel, assemble the spectra into a data cube, cor¬ 
rect for atmospheric dispersion, perform telluric 
corrections, and mosaic frames together. For some 
cases we utilized the Scaled Sky Subtraction mod¬ 
ule based on the technique outlined in Davied 
(120071 ). which scales the thermal continuum and 
OH line groups separately to provide optimal sky 
subtraction; we modified the module to include a 
smoother subt raction of the t herma l continuum, 
as described in lMedling et al. ( 2014 ). 


We imaged the tip-tilt stars for each galaxy 
during the course of the observations in order to 
obtain an estimate of the point-spread function 
(PSF). A Moffat function was fit to each tip-tilt 
star and then broadened according to the distance 
between the target and the tip-tilt star. This 
broadening accounts for the fact that the laser 
spot is not at infinity; therefore it only probes most 
of the turbulence that affects the target. As a re¬ 
sult, the further a target is from its tip-tilt star, 
the worse its PSF will be. The isokinetic angle for 
the Keck II AO system, at which the Strehl ratio is 
reduced to 37% of its peak value, is approximately 
75 arcseconds. We also use this an estimate of the 
Keck I AO system performance, which has not yet 
been characterized. 


1 Available at http://irlab.astro.ucla.edu/osiris/pipeline.html 
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Table 1 

Details of Observations 


Galaxy Name 

Redshift 

Pixel Scale 
(pc / 0.035") 

UT Date(s) 

Filter 

Exp Time on 
Target (minutes) 

Exp Time on 

Sky (minutes) 

CGCG436-030 

0.0315 

21.8 

2012 Sep 30 

Kn5 

30 

5 

IRASF01364-1042 

0.0493 

33.6 

2012 Oct 01 

Kbb 

20 

10 

IIIZw035 a 

0.0278 

19.5 

2011 Dec 10 

Kbb 

100 

50 

MCG+08-11-002 

0.0195 

13.9 

2012 Jan 02 

Kbb 

90 

50 

NGC 2623 

0.0196 

13.9 

2010 Mar 04, 2010 Mar 05 

Kn5 

100 

50 

UGC5101 a 

0.0413 

27.0 

2010 Mar 04, 2010 Mar 05 

Kn5 

80 

50 

Mrk273/UGC8696 b 

0.038 

26.4 

2012 May 22 

Hn4 

60 

30 

NGC 6240 N a 

0.0244 

17.2 

2009 Jun 17 

Kn5 

210 

75 

NGC 6240 S° 

0.0244 

17.2 

2007 Apr 21 

Kn5 

20 

10 

IRASF17207-0014“ 

0.0462 

29.8 

2011 May 23, 2011 May 24 

Hn4 

40 

20 


“Originally presented in IMedling et all 11201411 
b Originally presented in K! ct nil < 1201.4) ) 
c Originally presented in lMedling et all l i 201 ill 


3. Line Fitting 


4. Black Hole Masses 


We measured the kinematics following the 


metho ds presented in detail in Medline et al 
( 2011 . for stellar kinematics) and U et al. ( 20131 


fo r emission line kinem atics), and reviewed again 
in Medline et al.1 ( 2014h . We briefly describe these 


techniques here. 

For each measurement, we first calculate the 
signal-to-noise ratio in each pixel and b in th em us¬ 
ing op timal Voronoi tesselations (Cappellari & Copin 
120031 ) in order to require a certain signal-to- 
noise ratio for reliable measurements. We used 
thresholds of 20 and 3 per resolution element 
for stellar kinematics and emission line kinemat¬ 
ics, respectively. Once binned, stellar kinematics 
were fit using t h e pen alized pixel fitting routin(@ 
(|Cap pellari fe Emsellem l 12004!) and using A'-band 
templates of l ate-typ e giant s and supergiants from 
GNIRS ( Winge et al. 2009’). Emission lines were 
fit with Gaussian profiles to determine the flux, 
velocity, and velocity dispersion in each bin. For 
emission lines where multiple lines exist in a spe¬ 
cific band (the 5 H 2 transitions in the AT-band, or 
Br<5 and Bry), those lines were fit simultaneously, 
requiring that the velocity and velocity dispersion 
be consistent between lines. 


2 Available at http://www-astro.physics. ox. ac.uk/-mxc/idl/ 


Our black hole mass measurement techniques 
have been demonstr ated in two prior papers. In 
Medling et al.1 (|2011), we place limits on the mass 
of the black hole in the south nucleus of NGC 6240 
using 2-D stellar kinematic maps and two meth¬ 
ods. We calculate the lower limit to the mass by 
assuming the stars lie in a thin Keplerian disk, 
and fit a density profile including a black hole and 
a smoothly-varying spheroidal mass component: 
p(r) = Mbh + Po r ~ J - Because of the high spa¬ 
tial resolution of our data, we probe radii small 
enough that the first term dominates over the sec¬ 
ond. The measured black hole mass here is a lower 
limit because this method ignores velocity disper¬ 
sion (assuming that all measured velocity disper¬ 
sion is due to material along the line of sight, not 
intrinsic to the disk). Any intrinsic velocity dis¬ 
persion in the disk will be indicative of extra mass 
in the black hole. In the second scenario, we ob¬ 
tain an upper limit on the black hole mass by as¬ 
suming the opposite: that all measured velocity 
dispersion is intrinsic to the disk. Because the 
system is unrelaxed, there is unvirialized material 
along the line of sight increasing the velocity dis¬ 
persion but not knowing about the central black 
hole mass. To find this limit, we use Jean s Ax- 
isymmetric Mass models (JAM: [Cappellari]|2008 ). 
We use both methods here for all galaxies with 
measured stellar kinematics, except for the north 
black hole in NGC 6240. This nucleus contains 
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patchy dust that prevents a consistent light pro¬ 
file from being produced by the JAM models; in 
order to use this measurement technique on that 
black hole, equivalently high resolution imaging 
at a longer wavelength (L-band or beyond) is re¬ 
quired. 


In U_et ah (|2013l ). we measure the mass of 


the black hole in the north nucleus of Mrk273 
in a similar manner. Since no stellar kinematics 
were available, the JAM models were not feasible. 
When only gas kinematics are available, we model 
the gas as a thin Keplerian disk and proceed as 
in the first case above. The black hole mass mea¬ 
sured for this galaxy agrees well with the measure- 
rnent made using the_kinematics of an OH maser 
( Klockner fe Baanl l2004h . and therefore provides 
independent justification for this technique. We 
note that the thin Keplerian disk approximation 
provides reasonable Mbh agreement to the maser 
measurement e ven t hough the disks analyzed in 
Medline et al. ( 20141 ) are thick rather than thin, 


with v/<7 ~l-4. We note that this thickness may 
cause our black hole mass measurements to under¬ 
estimate the true black hole mass by up to ~ 25%. 


With this approach, we are measuring the un¬ 
resolved central mass within ~ 25 parsecs. This is 
likely dominated by the black hole, and we there¬ 
fore refer to it as such. However, this mass in¬ 
cludes the accretion disk of the black hole and 
may include a reservoir of gas feeding the accre¬ 
tion disk. If the accretion disk is fueled smoothly 
by the nuclear disk, this mass will be accounted 
for by the radial profile; we cannot account for a 
pile-up of gas below our resolution limit. The pos¬ 
sible implications for this will be discussed further 
in Section [5] 


In each case, errors on the black hole mass mea¬ 
surement were obtained using a Monte Carlo ap¬ 
proach. We created a model datacube and added 
random noise appropriate to the signal-to-noise of 
each spectrum and then refit the black hole mass 
for 100 iterations. The width of the resulting mass 
distribution is the error on the measurement. 


We have performed the above analysis for 9 ad¬ 
ditional black holes, listed in Table [I] and include 
the results in Table [2] and the fits in Appendix [A] 
As in previous papers, these methods produce 
kinematic models in reasonable agreement with 
observations for most cases. The ability to mea¬ 
sure black hole masses using different tracers, and 


additionally to use bracketing assumptions about 
velocity dispersion for stellar models, improves the 
robustness of these results. 


The numerical results from our black hole mass 
measurements are listed in Tabled along with pa¬ 
rameters relevant to black hole scaling relations, 
pulled from the literature: the stellar velocity dis¬ 
persion of the bulge, the stellar mass of the bulge, 
and the total stellar mass of the galaxy. Though 
consistent treatment of such parameters when 
drawing from different literature sources is diffi¬ 
cult, care was taken to select measurements most 
appropriate to the scaling relation. The bulge 
luminosities and total stellar masses were taken 
from the GOALS papers lHaan et~ah d 201 ll ) and 


U et al. ( 2012h . respectively. Haan et al. (|2011 ) 


measured the bul g e lum inosities using GALFIT 
( Peng et, al. I2002L l2010h to decompose iJ-band 
NICMOS images into multiple Sersic compo¬ 
nents. These galaxies are morphologically com¬ 
plex, which can produce an increased uncertainty 
in bulge fit parameters; th us, t he erro rs in bulge 
luminosities reported by Haan et al.1 ( 2011 1 and 
reproduced in Table [2] may be larger than the 
equivalent measurements from Hubble imaging 
of isolated galaxies. It is worth considering that 
we use the formal statistic al errors produced bv 
GALFIT; the analysis of Haussler et al.1 ( 20071 1 
has shown that these underestimate the true er¬ 
rors, which can contain contributions from profile 
mismatch and nearby neighbors, both releva nt to 
our sample. However, lHaussler et al. ( 2007h also 
find that GALFIT shows no systematic offset to 
fitting bulge luminosities of bright objects and 
that it is able to appropriately handle contami¬ 
nation by neigh bors when fit simultaneously (as 
was done in lHaan et all l201ll) . Still, for a more 
complete understanding of the errors associated 
with each bulge luminosity measurement, we re¬ 
fer the readers to t h e mo d els an d residuals for 
each system of lHaan et al.l ( 2011 1. published on¬ 
line in Figure Set 7. Our total stellar masses were 
measured from SEDs using photometry masks 
( U et al.1 2012t ) designed to incorporate all of the 
light from the galaxy. We are thus not underesti¬ 
mating the total stellar mass, as can happen with 
fixed-aperture photometry, and we avoid potential 
biase s related to varying ho st galaxy effective ra¬ 
dius (Hopkins et ahl 12007 afbL 20091 : Beifiori et al 


20121 b The GOALS survey papers allow us to 
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be confident that the bulge luminosities and total 
stellar masses were measured consistently across 
our sample. However, no such survey paper ex¬ 
ists currently to measure the stellar velocity dis¬ 
persions of the bulges. Instead, we selected the 
stellar velocity dispersion measurement from the 
aperture that most closely aligned with the size 
of the bulge. Unfortunately, no bulge information 
for MCG+08-11-002 is available in the literature. 
Although it cannot be placed on scaling relations 
at this time, we include our black hole mass mea¬ 
surements for completeness. 
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Table 2 

Measured Black Hole Masses 


Galaxy Name 

Tracer 3- 

Mbh 

^ ,bulge 

i / ^ H,bul qe \ 

lo s(— 

total 

Galaxy 



(M 0 ) 

(km s -1 ) 

(M 0 ) 

References 13 

CGCG436-030 

Br7 

4.59X 10 s 

175 ± 9 

10.84 ± 9.4 

5.5 X 10 9 

1,2,3 

IRASFO1364-1042 

Br7 

2.37+0 0 , 1 X 10 9 

103 ± 12 

10.42 ± 8.39 

5.6 X 10 10 

1,2,3 


h 2 

2.121“;°® x 10 9 





IIIZw035 

stars - disk 

> 6.8l°'J x 10 8 


11.15 ± 9.11 

1.8 X 10 10 

1,2 


stars - JAM 

< 2.0+°'® x 10 9 






h 2 

s.sgllo 0 -! 3 x 10 s 






Br7 

00 

O 

X 

■*< CO 

do 
+ 1 

cd 





MCG+08-11-002 

stars - disk 

> 8.7l l l X 10 7 






stars - JAM 

< 5.9 +7 'g x 10 8 






h 2 

6.9lJ'° X 10 7 






Br7 

2-3l l.l X 10 7 





NGC 2623 

stars - disk 

> X 10 8 

95 ± 13 

10.61 ± 8.57 

2.4 x 10 10 

1,2,4 


stars - JAM 

< 4.7j^;j? X 10 8 





UGC5101 

stars - disk 

> 6.5+21 X 10 8 

287 ± 11 

11.58 ± 9.75 

1.7 X 10 11 

1,2,5 


stars - JAM 

< 5.4t. 7 0 0 ' 4 8 x 10 8 





Mrk273 N 

[Fe II] 

1.0 ± 0.1 X 10 9 

285 ± 30 


6.9 X 10 11 

1,6 

NGC 6240 N° 

stars - disk 

> 8.8tg^ X 10 8 

174 ± 54 

10.81 ± 9.38 

3.9 X 10 11 

2,7,8 

NGC 6240 S 

stars - disk 

> 8.7 ± 0.3 X 10 8 

236 ± 24 

11.29 ± 10.25 

3.9 X 10 11 

2,7,8 


stars - JAM 

< 2.0 ± 0.2 x 10 9 





IRASF17207-0014 E 

[Fe II] 

2.5t° 0 0 3 3 x 109 

229 ± 15 

10.39 ± 9.06 

1.52 x 10 11 

2,6,7 

IRASF17207-0014 W 

[Fe II] 

8-9-1. 6 x 11)7 

229 ± 15 

10.39 ± 9.06 

1.52 x 10 11 

2,6,7 


a All measurements with gas tracers use the thin disk method, those with stellar tracers use the disk or JAM method as 
marked. 

b Galaxy parameters in co lumns 4-6 t a ken fr om t he following referen c es: 1 - IU et all 12013). 2 - IHaan et all J2011IL 3 
-IT^monti^^^^J|20^) , 4 - Shier et al.l 1199611 , 5 - IRothberg & Josephl 12U0f.il) , 6 - IPasvrE^^: al.l 12006h . 7 - iHoweh^^; al.l 


1201011. 8 -ITecza et aUllMom 

c Due to extensive and patchy dust coverage, a JAM model could not be completed for NGC 6240N without additional 
longer-wavelength (L-band or beyond) imaging data. 


























5. Mergers and Black Hole Growth 

5.1. Placing Black Holes on Scaling Rela¬ 
tions 

Using our measured black hole masses and host 
galaxy parameters from the literature, all listed in 
Table [21 we place the galaxies on several black hole 
scaling relations. 


5.1.1. The Mbh ~ Relation 


In Fi gure El we s how the Mbh — cr* rela¬ 
tion from [McConnell & Ma 820131 ) with our points 
overlaid. 



a, (km/s) 


Fig. 2.— The Mbh — <r* relation for isolated 
galaxies from McConnell fe Mai (l2013h (black data 
and best-fit line: log (Mbh/Mq) = 8.32 + 
5.64 log(cr*/200 km/s) with our merging galaxies 
marked in large colored points. Symbols indicate 
which method was used for that mass determi¬ 
nation: gas disk modeling in circles, stellar disk 
modeling in upward triangles, and JAM modeling 
in downward triangles. Stellar disk models rep¬ 
resent lower limits to the black hole mass, while 
JAM models represent upper limits; see text for 
details. Mrk273 N’s point (6; brown) is the mass 
measurement for only the north bla ck h o le, thoug h 
there is a second in the system fsee lU et al. 2013 ). 


The black holes in our sample he either within 
the scatter of the relation or above it. Although 
a particular galaxy may fall within the scatter of 


the relation, we want instead to consider whether 
our sample as a whole is behaving as predicted by 
the Mbh~&* relation. We calculate the distances 
from the Mbh — relation for both our dataset 
and the reference population, and then perform 
a two-sided Kolmogorov-Smirnoff test on the two 
populations. We find a p-value of 0.003, indicating 
that the chance of these two being drawn from the 
same population is less than 1%. To mitigate the 
possible effects of outliers in our small sample, we 
use the balanced bootstrap methocd to estimate 
the average offset from the relation. This reveals 
an average offset of 0.15 ± 0.06 in plot units; we 
thus consider this offset a 2.6-sigma result. 

Our sample of black holes are thus, on aver¬ 
age, above the relation, suggesting with moderate 
significance that black holes may grow in mergers 
before the bulges are virialized. Measuring the ve¬ 
locity dispersion of stars in the bulge can be diffi¬ 
cult to define when two bulges are in the process of 
merging. This variation in measured velocity dis¬ 
persion (both due to the process of merging and 
due to geometric variations based on the line-of- 
sight) has been quantified in the hy drod ynamical 
galaxy merger simulations of Sticklev fc Canalizol 
( 20141) . Their analysis shows a maximum pre¬ 
dicted mismeasurement in velocity dispersion of 
approximately 50%, which is insufficient to bring 
all our systems on to the relation. 


5.1.2. The Mbh — Lh, bulge Relation 

In Figure [3[ we show the Mbh — Lh, bulge rela¬ 
tion. The bac kground galaxies were pl otted using 
Lh, bulge from iMarconi fc Hunt ( 2003 ). with up¬ 
dated black hole masses from McConnell & Ma 
(|2013l ). We have updated the best-fit line with 


these numbers, producing the relation: 


to9< ^ ) = 8 - 22 + L 06 M 13 ^ ) 


The black holes from our sample fall either 
within the scatter or above the Mbh ~ Lh, bulge 
relation. As with the Mbh — 0 * relation, we per- 


3 Balanced bootstrap resampling is similar to ordi¬ 
nary bootstrap resampling but requires that each 
observed value appears with equa l freq uency in 
the resampled data (e.g. I Gleasonl 11988T) . We 

used the IDL routine BBOOTSTRAP, available at 

http://www.astro.Washington.edu/docs/idl/cgi-bin/getpro/library!4.html?I 
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form a two-sided Kolmogorov-Smirnoff test to de¬ 
termine if our sample as a whole falls significantly 
above the relation; we find a p-value of 0.018, indi¬ 
cating that these black hole masses are unlikely to 
be pulled from the sample distribution as the ref¬ 
erence sample. We also calculate the mean offset 
using the balanced bootstrap method, 0.47 ± 0.14 
in plot units, showing an offset with ~ 3.4-sigma 
significance. 


We note two possible calibration effects when 
considering our bulge luminosity measu rements. 
First, although iMarconi fc Hunt ( 20031 ) did not 
correct their bulge luminosities for internal ex¬ 
tinction, we consider the effects such a correc¬ 
tion would have on our sample because our galax¬ 
ies are quite dusty. Severa l of our galaxies hav e 
JHK imaging published by Scoville et al.l ( 2000h . 
enabling a direct estimate of typical extinction 
over the central few kiloparsecs. These calcu¬ 
lations estimate that the intrinsic luminosity of 
the bulges in our sample may be ~ 10% higher. 
This correction (indicated by red arrow in the top 
left of Figure El is smaller than the size of our 
points and can be safely ignored. Second, we in¬ 
vestigate whe ther the varying bu lge Sersic indices 
of the fits in lHaan et~ah ( 201lh may be biasing 
our luminosity estimates. When fixing n to 4, 
the bulge luminosities of several of our galaxies 
change significantly, some higher and some lower 
(S. Haan, private communication). The net ef¬ 
fect with this change is less scatter between our 
galaxies: a smaller but more statistically signif¬ 
icant offset. We conclude that, while the Sersic 
index fit may have a strong effect on an indi¬ 
vidual system’s predicted black hole mass when 
using the Mbh — Lh, bulge relation, the offset in 
our data is not caused by the fitting technique of 
these bulge luminosities. We also note that the 
two galaxies which fall furthest above the rela¬ 
tion, IRASF01364-1042 and IRASF17207-0014, do 
not show atypically large residuals in the GALFIT 
bulge fits (indeed, the residuals of the former ap¬ 
pear smaller than average for this sample), which 
suggests that their of fset is not du e to er rors in the 
bulge luminosities of Flaan et al . fall ill . 


Again, our sample of merging galaxies falls 
above this relation, suggesting that either stars 
have not been forming in the bulge as quickly as 
the black holes have been growing, or that a large 
mass of stars is still strewn about in tidal tails and 


unrelaxed features that will eventually become 
part of the bulge. As the H -band traces older pop¬ 
ulations of stars, it may be that some stars have 
recently formed in the bulge, but that this mea¬ 
surement misses them. Although Mbh — Lbuige 
relations do exist for bluer colors, we are unfortu¬ 
nately unable to utilize them for this sample be¬ 
cause the large amounts of dust obscure the galaxy 
cores, making bulge luminosity measurements in 
the optical highly uncertain. We have compared 
our bulge luminosities to those fit by Kim et ah 
( 2013 ) in the /-band, and found that (for typical 
mass-to-light ratios) our //-band bulges are more 
massive. This simply confirms that dust is a big¬ 
ger problem even in the /-band than missing stel¬ 
lar populations are in the //-band for this sample. 


5.1.3. The Mbh — M* Relation 


In Figure IH we p l ot Mb h — M* using points 


from Bennert et al.1 ( 20111) and ICisternas et al 


(|2011al ). These authors suggest that higher red- 
shift galaxies sit on the local Mbh — M*^uige 
relation if you include total stellar mass instead 
of only bulge mass. In this way, they attribute 
evolution in the black hole scaling relation to the 
evolution of the fraction of stellar mass in the 
bulge. We include as a solid line the upd ated 


fit to the Mbh — M*^ulge from McConnell & Ma 


( 2013 ). and plot our black hole masses in the same 
colored points as in Figures [2] and [3] 


To determine if this sample of systems is in 
agreement with the Mbh — M* relation, we per¬ 
form a similar two-sided Kolmogorov-Smirnoff sta¬ 
tistical test as above. We find th at, com p ared 
to the local elliptical galaxies of lMcConnell fc Ma 
(l2013h . the p-value obtained is only 0.01, less 
significant than above. When comparing to the 
higher-redsliift samples of total stellar mass, we 
obtain a p-value of 0.0008, and to the combined 
sample, a p-value of 0.004. The balanced boot¬ 
strap method finds an offset from the McConnell 
& Ma relation of 0.46 ± 0.14 in plot units, similar 
to the Mbh — Lbuige relation. 

We thus conclude that our sample of galaxies 
on average lies above this scaling relation as well. 
We also note that the one point that falls below 
the relation, the brown circle Mrk273 N (labeled 
6), represents the mass in only one out of two 
black holes and is therefore artificially low. Be¬ 
cause this relation takes into account the entirety 
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of stellar mass, the fact that many of our points 
fall above the relation suggests that star forma¬ 
tion as a whole is not outpacing black hole growth 
during these mergers. 


5.1.4- Overall Significance of Scaling Relation 
Offset 


In the previous sections, we have compared 
our galaxies with the Mbh — <J* and Mbh — 
M* t bul ge relations from McConnell fc Ma ( 2013 '). 
We have also compared them to those from 
Kormendv fc HcJ ( 2015 '). and find that our offset 


remains. 


Two galaxies appear to drive the offset in 
the Mbh — c* relation: IRASF01364-1042 and 
NGC 2623. Although we note that these are not 
outliers in a statistical sense, we must be certain 
that they are not unfairly biasing our results. Be¬ 
cause our sample for each scaling relation has only 
8 points, throwing out these two galaxies has a 
significant effect in all cases through reducing the 
sample size. We note, however: when removing 
the two galaxies from the Mbh — cr* relation, the 
measured offset decreases quite a bit, while the 
uncertainty in the offset also increases (producing 
only a 1.9-er result). For the other two relations, 
the measured offset decreases only slightly, and 
the main effect is an increase in uncertainty (so 
the significance still decreases from 3.3-cr to 2.3-cr). 
This demonstrates that the offsets aren’t driven 
by these two galaxies in the Mbh — Lbuige and 
Mbh ~ M* relations. 

Due to the small sizes of our samples, the ap¬ 
propriate way to deal with potential outliers is to 
use the balanced bootstrap test, which we have 
done. The balanced bootstrap performs hundreds 
of resamplings in order to estimate the significance 
of a result, giving equal weight to all measure¬ 
ments. (Non-‘balanced’ bootstrapping may re¬ 
sample some points more than others, which could 
cause a bias when one or more potential outliers 
exist. By balancing the bootstrapping, we force 
the equal treatment of all measurements; this pre¬ 
vents a single point (or two) from biasing the re¬ 
sults.) By doing so, we are able to handle outliers 
(or near-outliers) without diminishing our sample 
size. 


The black holes in this sample of gas-rich ma¬ 
jor mergers consistently fall on or above the scal¬ 


ing relations for passive galaxies with stellar ve¬ 
locity dispersion, bulge luminosity and total stel¬ 
lar mass. Although a number of our galaxies lie 
within the scatter of these scaling relations, our 
statistical tests show they are significantly more 
likely to fall on the upper side of the scatter than 
randomly distributed about the relations. Though 
the offset from the Mbh — c* relation is only a 2.6- 
sigma result, the agreement with the two other 
significant relations lends it confidence. Taken 
together, these results show that increased black 
hole growth in these mergers has already begun 
and may be proceeding more rapidly than growth 
of the host galaxy. 

Although we have a relatively small sample of 
black holes, it is striking that we already see a 
significant offset. As more black hole masses are 
measured in gas-rich mergers, the larger sample 
will provide a more accurate picture of how offset 
such systems are from black hole scaling relations. 


We caution the reader against considering the 
distance from a scaling relation as a direct proxy 
for merger stage. This is clearly evident when 
comparing where a system falls on one scaling re¬ 
lation versus another. For example, IRASF17207- 
0014 falls closer to the Mbh — u* relation than 
to the Mbh — Lh, bulge relation; for NGC 6240S, 
the reverse is true. This discrepancy could be 
caused by at least two possible mechanisms: 1) 
The host galaxy parameters trace different aspects 
of galaxy growth, which may evolve differently. 
For example, the bulge velocity dispersion traces a 
convolution of mass growth with the virialization 
timescale whereas the bulge luminosity evolution 
may rely less stringently on the dynamical relax¬ 
ation ti me. 2) As discu s sed ab ove, the simulations 
of lSticklev fc Canalizol ( 2014h show that the mea¬ 
sured stellar velocity dispersion of similar mergers 
oscillates during the course of the merger. Analo¬ 
gous simulations would have to be performed for 
bulge luminosity and total stellar mass measure¬ 
ments to determine if they ought to oscillate as 
well. Even if all three quantities were expected to 
oscillate during a merger, it seems contrived to re¬ 
quire that they oscillate in lockstep. To determine 
precisely how galaxies evolve along a scaling rela¬ 
tion (and if they evolve the same way along each 
scaling relation), a larger sample of systems is re¬ 
quired. By binning over merger stage, the range of 
distances from scaling relations can be seen, and 
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key parameters affecting this can be identified. We 
leave this analysis for future studies. 

5.2. Interpreting Overmassive Black Holes 

In this section, we discuss several plausibility 
arguments related to the above results. We con¬ 
sider the growth timescales and the available fuel 
for both black holes and the host galaxies. 


5.2.1. Has There Been Enough Time for these 
Black Holes to Grow? 


We posit that these systems followed black hole 
scaling relations before the merger began and that, 
since then, mass has been funneled to the center 
to enable black hole growth. For simplicity, we 
assume that each galaxy previously fell precisely 
on every scaling relation; of course, as each re¬ 
lation contains intrinsic scatter, this is unlikely. 
However, we have no reason to believe that these 
systems were systematically offset from scaling re¬ 
lations before the mergers began, and thus this 
simplification should average out in the order-of- 
magnitude calculations of this and the following 
sections. 

With that assumption, we consider how long it 
would take for a black hole, growing at the Ed¬ 
dington rate, to reach this far beyond its original 
mass. In Table [3j we calculate this for our sam¬ 
ple of galaxies. We see that the required times for 
these black holes to grow, while accreting at the 
Eddington rate, are a few tens to hundreds of mil¬ 
lions of years. As this is less than or comparable to 
a merger timescale, it is reasonable that these sys¬ 
tems could have lain on scaling relations before the 
merger, and risen only after the merger began. As 
black hole growth is known to be episodic, we note 
that these systems could often have been accret¬ 
ing at sub-Eddington rates for timescales less than 
those noted, and that the accretion times need not 
have been continuous. 


The merger of one o f these object s . NGC 2623, 
has been modeled by Privon et al.l ( 2013 ') using 
H I kinematics from the Very Large Array. Their 
dynamical models indicate that this object is 2.2 x 
10 s years past first pericenter passage, within the 
range of growth timescales we see for this object 
in Tabled (0.6 — 2.3 x 10 8 years). 

We also repeat that our black hole mass mea¬ 
surements produce dynamical masses of a central 


point source, which to our spatial resolution could 
be anything smaller than ~ 25pc. It is therefore 
plausible that some of the mass measured is in 
the accretion disk (or, perhaps, in a nearby reser¬ 
voir feeding the accretion disk). If that’s the case, 
any mass that has not accreted directly onto the 
black hole would not cause the black hole to turn 
on as an AGN. In such a scenario, the Eddington 
rate may not be a relevant limiting factor in this 
growth phase. 


5.2.2. How Long Would It Take To Return to 
Scaling Relations? 


The gas-rich major mergers we’ve studied here 
appear to have black holes significantly larger than 
would be expected based on their host galaxy 
properties. We assume that these systems, post¬ 
merger, ought to fall on scaling relations. As our 
systems are currently undergoing enhanced star 
formation, it is likely that the host galaxies will 
indeed continue to grow over the course of the 
merger, moving back towards the local scaling re¬ 
lations in Figures[2]m In Tabl effl we show the rele- 
vant star formation rates from Howell et al . H 20 T 0 I) 
and calculate how long it would take at this star 
formation rate to bring the total stellar masses 
up to those predicted by the measured black hole 
masses. 

We see the typical timescale for these ob¬ 
jects to reach predicted total stellar masses is 
1-2 Gyr at the present star f or matio n rate . The 
dynamical models of Privon et al. ( 20131 ) find 
merger timescales for their sample of 4 systems 
to range from 140 Myr - 1.4 Gyr, depending on 
the initial conditions of the merger. Apart from 
IRASF10364-1042, our galaxies are not far from 
this range. Though star formation is a variable 
quantity, it is plausible for these galaxies to return 
to the Mbh ~ Af* relation by the time they have 
finished merging. 


As noted above, it is possible that not all the 
mass we measured has been accreted by the black 
hole yet. If this is the case, it allows for the pos¬ 
sibility that some of this mass might be ejected 
via AGN-driven winds and may indeed never be 
accreted; we discuss this scenario in Section 15.31 
If this were to happen, the stellar mass growth 
need not be so high to return these systems to the 
M B h ~ M* jotai relation. 
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Table 3 

Black Hole Growth Timescales 


Galaxy Name 

Mean 

Mbh (Mq) 

From cr* 

, bulge 

From Lb, 

Implied Mass 
Growth (Mq) 

bulge 

Time at 
Edd(yr) 

From M „ 

, total 

Implied Mass 
Growth (Mq) 

Time at 
Edd (yr) 

Implied Mass 
Growth (Mq) 

Time at 

Edd (yr) 

CGCG436-030 

4.6e8 

3.6e8 

7.7e7 

2.8e8 

4.6e7 

4.5e8 

1.8e8 

IRASF01364-1042 

2.3e9 

2.3e9 

3.1e8 

2.2e9 

1.8e8 

2.2e9 

1.3e8 

IIIZw035 

2.6e8 



- 

- 

2.1e8 

8.5e7 

NGC 2623 

3.0e8 

2.9e8 

2.3e8 

1.9e8 

5.3e7 

2.3e8 

7.6e7 

UGC5101 

5.5e8 

- 

- 

- 

- 

4.9e7 

4.6e6 

Mrk273 N 

1.0e9 

- 

- 



- 

- 

NGC 6240N 

8.8e8 

7.8e8 

I.le8 

7.1e8 

8.2e7 



NGC 6240S 

8.9e8 

3.6e8 

2.6e7 

3.5e8 

2.5e7 

o.oeo 


IRASF17207-0014 

2.6e9 

2.1e9 

8.8e7 

2.5e9 

1.9e8 

2.1e9 

8.8e7 


Note.—I mplied mass growth columns give the difference between the measured Mbh and that predicted from the indicated 
host galaxy properties and scaling relations. The time at Eddington rate columns are caculated from Mb h, now / Mbh, predicted, = 
e t / T Sai ? where rs a l is the Salpeter timescale of 5.7 X 10 7 years. Columns with blank entries do not have measured host galaxy 
properties; entries marked with a — indicate systems that lie on or below scaling relations, and therefore require zero growth. 


Table 4 

Star Formation Requirements 


Galaxy Name 

AM., W / 

SFR (M 0 /yr) b 

Time Required (yr) c 

Available Mh 2 (-^©) d 

CGCG436-030 

1.5ell 

85.87 

1.7e9 


IRASF01364-1042 

6.7ell 

122.61 

5.5e9 


IIIZw035 

7.3el0 




NGC 2623 

7.9el0 

69.19 

l.le9 

1.66e9 

UGC5101 

1.6el0 

180.18 

8.7e7 

2.65e9 

Mrk273 N e 

- 

- 

- 

- 

NGC 6240 N+S 

1.7ell 

148.44 

1.2e9 

1.18el0 

IRASF17207-0014 

6.6ell 

501.22 

1.3e9 

7.13e9 


a A M* t totai gives the difference between the current M*,total (column 6 in Table m and that 
predicted from the current mean measured Mbh and the relevant scaling relation. 

b Star formation rates taken from [Howell et al.l l|2010f) . 

c Time required = A M* f total / SFR 

d From I Wilson et al.| J2008D when available 

e Mrk273 N already falls on the scaling relation, and therefore no additional star formation is needed. 
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5.2.3. Is There Enough Gas Mass Left to Grow 
the Host Galaxy? 

Regardless of the rate of star formation dis¬ 
cussed in Section 15.2.21 we might ask if there is 
enough molecular gas present in these galaxies 
to increase the sufficiently if all the gas 

were to be turned into stars eventually. We in¬ 
clude the available mass in cold gas measured for 
each of these galaxie s, when avail a ble, u sing H 2 
mass es measured bv IWilson et al.1 (1200811 . How¬ 
ever, Wilson et al.1 (2008) point out that the CO- 


to-H 2 conversion factor used is a primary source of 
uncertainty, and caution that their masses could 
be underestimated by up to a factor of 5, and that 
AGN contamination in UGC5101 could cause fur¬ 
ther underestimation. Taking these uncertainties 
into account, the gas masses could be sufficient 
in UGC5101 and NGC 6240, but perhaps not in 
NGC 2623 or IRASF17207-0014. 

We note that this estimate ignores atomic gas 
which may eventually cool and form molecular gas 
and then stars; this quantity could be substantial, 
and that these systems could also continue to ac¬ 
crete gas from the circumgalactic medium. How¬ 
ever, a number of ULIRGs are also hosts to mas- 


sive molecular outflows (IRupke fe Veilleuxl 2011 


Sturm et al.l2011 : Spoon et al.l2013tlVeilleux et al 

2013t h which could instead further deplete the gas. 

Thus, the low gas mass estimates may support 
the possibility that not all excess mass measured 
from the central regions will eventually be accreted 
to the central black hole. 

5.3. Are These Masses Truly Measuring 
Black Holes? 

Our technique produces a dynamical measure¬ 
ment of the mass produced by a point source in 
the centers of these galaxies. Though the most 
straightforward interpretation of this is that the 
mass measured describes a black hole mass, it is 
possible that the mass distribution merely appears 
as a point source to our spatial resolutions. That 
is, anything significantly smaller than our pixel 
scales of ~ 25pc, and above the disk’s radial den¬ 
sity profile, would appear as a black hole to our 
measurements. 


rule out this interpretation for two reasons: 1) nu- 
clear star clus ters commonly have radii ~ lOOpc 
( Murravl2009f) . which would not appear as a point 
source to our dynamical modeling, and 2) the 
continuum luminosity present is not sufficient to 
account for a significant fraction of the central 
masses being due to stars. (Full spectral en¬ 
ergy distribution fitting to the star clusters in 
and around the nuclear regions of these galaxies 
will be presented in a future paper.) We cannot, 
however, rule out the idea that gas present in the 
center could someday form a nuclear star cluster. 

Galaxy merger sim ulations commonly assume 
Bond i -Hoyle accretion dBondi fc Hovlell944 ; Bondi 
1952t ISpringel et al.1 2005bl) . a prescription in 
which effectively gas particles which enter within 
the sphere of influence of the black hole are auto¬ 
matically accreted so long as the Eddington rate 
is not exceeded. However, gas inflow likely has an¬ 
gular momentum which must be discarded before 
accretion can occur, which can delay the accretio n 


event by e.g. the viscous ti mesca le ((Power et al 
2011 : Wurster fc Thacker 2013bl l. Though we 


don’t have strong constraints on the properties 
of AGN accretion disks, estimates predict that 
viscosity could slow down accretion by more than 
a Hubble time (King 2008 1. Although the exis¬ 


tence of AGN confirms that accretion can occur 
on shorter timescales, the delay may not be negli¬ 
gible. This allows for the possibility that gas can 
build up in the core of a galaxy around the black 
hole without having yet been accreted. This ac¬ 
cretion disk (and, possibly, a reservoir feeding the 
accretion disk) of gas could contain a substantial 
amount of mass en route to the black hole. This 
mass would all appear as a point source to our 
dynamical models. 

Given the scaling arguments in the sections 
above, we also note the possibility that this pu¬ 
tative accreting gas may never make it to the 
black hole. Although the black hole dominates 
the gravitational pull on this gas, we find it plau¬ 
sible that when the black hole begins to accrete, 
its feedback may blow some of the rest of the 
gas away. Thus, the extra mass that we see 
here could be the reservoir for massive outflows 


Alatalo et al 


As discussed in Medline et al. 

20111 and 

2011: 

Ruoke & Veilleux 2011 

Sturm et al. 2011; 

U et al. 

20131), it is possible that some mass mav 

Veillcux et al.l 2013| Spoon et al. 2013). At mass 


be due to a nuclear star cluster. However, we 


outflow rates of tens or hundreds of solar masses 
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per year, the excess mass could easily be ejected 
in the time remaining in the merger. It is not 
currently known where the mass carried in these 
outflows is drawn from, or in particular, how much 
of it might come from within ~ 25 pc. 


Simulations have also suggested that nuclear 
gas disks may leave behind remnant stellar disks 
on 1-1 0 parsec sc ales with masses 0.1-1 Mbh 
( Hopkins fc Quataert 2010bl) ; if this were to oc¬ 
cur, the final black hole masses in our sample will 
be decreased by 10-50%. We note, however, that 
the black hole masses from which the scaling rela¬ 
tions were calibrated did not have sufficient resolu¬ 
tion to separate out remnant disks on 1-10 parsec 
scales either; thus we only introduce a bias if our 
sample is more likely to host such disks than the 
typical galaxy. We note that such remnant disks 
were discovered in galaxy merger simulations but 
that they did not form until after the final coa¬ 
lescence of the two nuclei. Thus it is not clear 
whether we would expect such disks (broken from 
the larger nuclear disk on hundred parsec scales) 
to be present in our sample. 


Though the above possibilities may moderate 
the black hole masses measured by attributing 
some mass to other features, we still conclude 
that the black holes (combined with their accre¬ 
tion disks) in this sample of gas-rich mergers do 
not lie below black hole scaling relations, and in 
many cases are significantly more massive than 
their host galaxy properties would predict. This is 
directly opposite the predictions of quasar-mode 
feedback theories, in which black hole growth is 
delayed until the final stages of the galaxy merger 
when it cuts off star formation. As we note that 
(U)LIRGs are prototypically expected to evolve 
into quasars ( Sanders fc Mirab el 19961 1. we find 
this as further evidence of some delay in feed¬ 
back timescale during a merger. Such a delay 
(due e.g. to a viscous accretion disk, as discussed 
above) could reconcile the late-stage AGN feed¬ 
back paradigm with the early mass growth pre¬ 
sented here. Testing this scenario requires spatial 
resolutions higher by a factor of 10-20, which may 
be feasible in the upcoming era of thirty-meter- 
class telescopes. 


5.4. Comparison with Previous Findings 


These results differ from the conclusions drawn 
by Kormenclv & Hoi (2 0131 ), who show five merging 


galaxies falling below the Mbh — cr* and Mbh — 
Lk, bulge relations. They suggest that black hole 
growth might lag bulge formation in mergers, the 
opposite of the findings of this work. It is impor¬ 
tant to note the differing samples when consider¬ 
ing these discrepant results. Our sample of merg¬ 
ing galaxies contains major gas- rich mergers; the 
mergi ng galaxies considered in iKormendv fc Ho 
( 2013 ) are mainly minor mergers and/or gas-poor 
mergers. It is not unreasonable to imagine that 
such different types of mergers would have differ¬ 
ent effects on their central black holes. For ex¬ 
ample, a minor or gas-poor merger might use up 
all incoming new gas in star formation before it 
reaches the central black hole, causing a delay in 
black hole accretion until stellar feedback produces 
new gas to fuel it. A major gas-rich merger may 
have sufficient gas inflow to enable both star for¬ 
mation and black hole growth to begin relatively 
early. 


A sample of black holes in submillimeter 
galaxies appear to fall below the locally-defined 
Mbh — M*.huinp. relation by 1-2 orders of mag - 
nitude ( Borvset al.l 120051 : 1 Alexander et al.l 12008 1. 
These galaxies are thought to be the high-redshift 
(z~ 2) counterparts of local ULIRGs and therefore 
might be reasonable to expect that they would be¬ 
have in the same way as the sample presented here. 
However, the discrepancy may in fact be due to 
their usage of broad line velocity dispersions to in¬ 
fer Mbh ; indeed, wh en applying the same me thod 
in one local ULIRG, Alexander et al. ( 2008 ) find 
rough agreement with their submillimeter galax¬ 
ies. It is also possible that differences between 
the samples such as Eddington ratio, redshift, and 
merger stage might affect the black hole growth 
timescales. 
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Fig. 3.— The M BH - L H}bu i ge relation. 

The black point s have luminosities taken from 


Marconi fc Hunt! ( 2003) and updated black hole 
masses from McConnell fc Mai ( 2013h . The 
solid line is a log-log fit to the updated data: 

) = 8-22 + 1.06 l 0 g(±jtfs%). Galax¬ 
ies presented here are indicated by large colored 
points, with parameters listed in Table [2] The 
bulge luminosities were not corrected for extinc¬ 
tion; the arrow in the top left shows the typical 
magnitude of such a correction (10%). Symbols 
indicate which method was used for that mass de¬ 
termination: gas disk modeling in circles, stellar 
disk modeling in upward triangles, and JAM mod¬ 
eling in downward triangles. Stellar disk models 
represent lower limits to the black hole mass, while 
JAM models represent upper limits; see text for 
details. 
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Fig. 4.— The .Mb#—M* relation. The grey points 
have luminosi t ies an d t otal stellar masses from 
Bennert et al. ( 201ll) and lCisternas et all ( 2011 a I. 


who suggest that higher redsliift galaxies sit on 
the local M BB — M*,bulge relation if you include 
total stellar mass instead of only bulge mass. In 
this way, they attribute evolution in the black hole 
scaling relation to the evolution of the bulge/disk 
fraction with redshift. We include also black 
points for local elliptical gala xies an d the fit to 
for Mbh ~ bulge from ( McConnell fc M«i1:2013 . 
log (M B h/Mq) = 8.46 + 1.05 log(M*/lO n M 0 ). 
Merging galaxies from this work are indicated by 
large colored points, with parameters listed in Ta¬ 
ble [2l Symbols indicate which method was used 
for that mass determination: gas disk modeling in 
circles, stellar disk modeling in upward triangles, 
and JAM modeling in downward triangles. Stel¬ 
lar disk models represent lower limits to the black 
hole mass, while JAM models represent upper lim¬ 
its; see text for details. NGC 6240’s two nuclei 
were summed together and appear as one point 
(7; dark purple) on this plot, paired with the en¬ 
tire system’s total stellar mass. Mrk273 N’s point 
(6; brown) is the mass measurement for only one 
black hole, though there is a second in the system 
(see ITT et a,1.11201 3). 
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6. Conclusions 


We have presented high spatial resolution near- 
infrared integral field spectroscopy of nuclear disks 
in nearby merging (U)LIRGs. These gas-rich 
mergers funnel gas to the galaxy centers where 
the gas forms disks; stars then form in situ, cre¬ 
ating a stellar disk as well. T hese dis ks, which 


were described in detail in iMedling et al.l (120141 ! , 


provide an alternative method to measuring black 
hole masses. Three-integral orbital superposition 
models rely on an assumption of virialization and 
emission line black hole mass diagnostics require 
broad lines; neither of these methods is feasible in 
such systems. However, with adaptive optics, we 
have produced two-dimensional kinematic maps at 
spatial resolutions close to or better than radius of 
the sphere of influence of the black hole. Because 
the dynamical timescale this close to a black hole 
is short, the material quickly becomes virialized, 
and kinematic modeling is possible. 

We implemented two different techniques us- 
ing these kinematic maps , origina l ly pre sented in 
Medling et ah | (l201llf and lll et al.1 (I2013I! . to mea¬ 


sure 9 new black hole masses, producing a sample 
of black hole masses in 9 merging systems. By ob¬ 
taining both gas and stellar kinematics for many 
of these systems, we are able to provide multiple 
independent measurements for a large fraction of 
our sample. Gaseous disks are modeled as thin 
Keplerian disks. Stellar disks are more likely to 
have pressure support, but measuring this is not 
straightforward, so we place bracketing measure¬ 
ments using two assumptions. To provide a lower 
limit to the black hole mass, we assume the stars 
are also in a thin Keplerian disk; any pressure sup¬ 
port would increase the black hole mass. We also 
use Jeans Axisymmetric Mass models to take into 
account the velocity dispersion; since this velocity 
dispersion may be inflated by intervening unviri- 
alized material along the line of sight, this mass 
measurement is an upper limit. 


As the black hole masses measured here are 
in fact masses enclosed in ~ 25pc, some of the 
mass attributed to black holes here may be gas 
that has not yet accreted onto the black hole. Al¬ 
though this mass is gravitationally bound to the 
black holes, we note that it is possible that future 
AGN-driven outflows may carry a fraction of it 
away, thus diminishing the final black hole mass. 


Though such outflows have now been seen in a 
number of (U)LIRGs, further studies will reveal 
how frequently they occur, how much mass they 
carry, and how much of the loaded mass comes 
from within ~ 25pc. 


We have placed these enclosed masses on three 
local black hole scaling relations, Mbh — o - *, 
Mbh-Lh, bulge, and using host galaxy 

parameters from the literature. We find that 
though individual galaxies often fall within the 
scatter, as a population they are significantly off¬ 
set towards higher black hole masses than isolated 
systems. Several systems have central masses that 
sit considerably above scaling relations. 

Under the assumption that these galaxies fol¬ 
low black hole scaling relations before and af¬ 
ter their mergers, this suggests that black holes 
grow more quickly than their host galaxies dur¬ 
ing a major gas-rich merger. This does not line 
up with theoretical reasoning that star fo rma¬ 
tion h as easier acces s to inflowing gas (e.g. ICen 
20121: iHonkinsI l2012h or that a final burst of 


black hole growth is likely to shu t off star for¬ 


mation (“quasar-mode fee dback”; iHoukins et al 


2005 : Springel et al.1 [200 5a!. unless that feedback 


i s delayed due to viscosity in the accretion d isk 
( Power et al. 2011 : Wurster fe Thackeill2013bll or 
similar. We conclude that black hole fueling be¬ 
gins early in a gas-rich merger, and can initially 
outpace any simultaneous starbursting or bulge 
formation. 


These results have focused on a sample of 
gas-rich mergers in late stages of merging. To 
trace black hole growth through the entire merger 
sequence, a larger sample of black hole masses 
should be obtained for systems at a range of 
merger stages. Such an observational sample 
could provide the basis for distinguishing between 
prescriptions for black hole accretion and feed¬ 
back^ physics_in ga laxy m erger simulations (e.g. 


Wurs ter fc Thacker 2013a!. 
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A. Dynamical Mass Models 
A.l. Fits and Residuals 


In this section, we show the best-fit dynamical mass models for each galaxy for each modeling technique 
in Figures 1511221 As discussed in Section[4j each black hole was dynamically fit using a subset of the following 
methods: thin disk fitting to Bry, Pa a, or [Fe u] velocity map, thin disk fitting to H 2 velocity map, thin 
disk fitting to stellar velocity map, JAM modeling to stellar velocity and velocity dispersion maps. Galaxies 
with multiple tracers show reasonable agreement in measured masses, which can be seen in Table [2] and 
Figures [2] [U Residuals from each fit, divided by the measurement error, are included in these figures and 
demonstrate the goodness-of-fit of each model. The fitted disk parameters are also given in Table [5] 


We note that many galaxies have complex kinematic structures that deviate from pure rotation in one or 
more of the tracers. These can be due to tidal streams or spiral arms along the line-of-sight, or to inflow 
or outflow when seen in gas tracers. Where possible, we have masked these out to avoid biasing the disk 
models. Masked spaxels appear white in the figures. A full analysis of the deviations from pure rotation and 
subsequent implications for inflow or outflow will be presented in a future paper. 

The positions of the black holes, plotted at (0,0) in each of the following figures, were determined primarily 
from the kinematic centers rather than the photometric centers. However, in the case of IRASF17207-0014, 
the kinematics are sufficiently compl ex that the two black holes were fixed to the photometric center of each 
disk measured with GALFIT in lMedling et al. ( 20l4) . 
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Fig. 5.— a) Map of Br 7 velocities of the inner region of CGCG436-030. b) Map of model velocities of best 
fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and b), 
velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 6.— a) Map of Paa velocities of the inner region of IRASF01364-1042. b) Map of model velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 7.— a) Map of H 2 velocities of the inner region of IRASF01364-1042. b) Map of model velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 8.— a) Map of stellar velocities of the inner region of IIIZw035. b) Map of model stellar velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 9.— Left: Map of the symmetrized v rms = (v 2 + a 2 ) 0 5 of IIIZw035. Center: Model v rms of galaxy with 
black hole. Right: Residual map of model - data divided by errors in v rms . All three panels are centered on 
the black hole and have been rotated so the major axis is horizontal. 
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Fig. 10.— a) Map of Bi '7 velocities of the inner region of IIIZw035. b) Map of model velocities of best 
fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and b), 
velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 


A.2. Black Hole Spheres of Influence 

One key feature of our thin disk black hole mass models is the inclusion of a radially-varying mass 
component to incorporate the mass profile of the nuclear disk. (Our JAM models also incorporate galactic 
mass, but by fitting a mass-to-light ratio.) We note that our fits use only a single power-law fit to the mass 
profile, so any pile-up of matter below our resolution limit will be considered part of the point mass. We 
include the mass profiles measured here for completeness, and for the purposes of calculating the implied 
radius of the sphere of influence of each black hole. Recall that the mass profile is fit with the form 
M(r) = por 1 . In Table Owe report the normalization and power law fit for each of these parameters. For 
clarity, instead of reporting poi we report the mass in the disk enclosed at 100 parsecs. 

A black hole’s sphere of influence is defined as the region within which the black hole contributes 50% of 
the mass. We can therefore use our thin disk mass profiles to calculate the radius of the implied sphere of 
influence. This is interesting both for the purposes of understanding black hole dynamics and for confirming 
that we indeed have sufficient resolution to measure the black hole masses appropriately. 

As the fitted mass profiles of the disk and black hole masses are not independent, there is an issue of 
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Fig. 11.— a) Map of H 2 velocities of the inner region of IIIZw035. b) Map of model velocities of best fit 
black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and b), 
velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 


covariance to consider. That is, when a model fits a higher black hole mass, it will likely put less mass in the 
disk, and vice versa. When calculating errors on each of these parameters, we marginalize over the others: 
our black hole mass errors are calculated from the width of the black hole mass distribution in our Monte 
Carlo simulation, taking into account all possible fitted disk profiles. Similarly, we report our errors on disk 
profile parameters using the widths of the distributions taking into account all fitted black hole masses. This 
is an appropriate way to handle covariance when drawing conclusions about the black hole masses, as is done 
in this paper. To account for the covariance issue in the measurement of ?y n /z, we calculate ri„fi for every 
fit in the Monte Carlo simulation and draw our conclusions based on that. 

In one case (IRASF01364-1042), the black hole mass is comparable or larger than the entire nuclear disk. 
The calculated radius of influence is thus artificially large, because presumably the mass profile of the galaxy 
changes beyond our field of view. The H 2 fit to IIIZw035 is also artificially large; we attribute this to 
streaming motion in the H 2 emission compromising the larger scale fit. This and other evidence of outflows 
will be presented in an upcoming paper. The affected radii are reported in parentheses in Table [5] 

We see that, for our sample, we clearly resolve the spheres of influence of these black holes. We note that 
this method of calculating the size of the sphere of influence is uncommon, as a full mass profile is required 
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to determine where M{r) = Mbh■ Our capacity to resolve these scales exceeded our expectations because 
the black holes are more massive than predicted. 
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Fig. 12.— a) Map of stellar velocities of the inner region of MCG+08-11-002. b) Map of model velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 13.— Left: Map of the symmetrized v rms = ( v 2 + a 2 ) 0 - 5 of MCG+08-11-002. Center: Model v rms of 
galaxy with black hole. Right: Residual map of model - data divided by errors in v rms . All three panels are 
centered on the black hole and have been rotated so the major axis is horizontal. 


Table 5 

Measured Disk Parameters 


Galaxy Name 

Tracer 

Mbh 

(Mq) 

Mgal(< lOOpc) 
(M 0 ) 

7 

rinfl 

(pc) 

rinfl 

(pixels) 

CGCG436-030 

Br7 

4.59±g'51 X 108 

6.0 ± 3.5 X 

10® 

0.7 ± 0.2 

390 ± 60 

17.9 ± 2.8 

IRASF01364-1042 

Br7 

2.37t° ( j 0 1 1 X 10 9 

6.2 ± 1.1 x 

10 4 

2.0 a 

(2.3 ± 0.1 x 10®) b 

(7.0 ± 0.3 x 10 4 ) 


H 2 

2.12±° °l x 10 9 

1.0 ± .1 X 

10® 

2.0 

(2.0 ± 0.1 x 10®) 

(6.0 ± 0.3 x 10 4 ) 

IIIZw035 

stars - disk 

> 6.8i 9 J X 10 s 

2.7 ± 0.4 x 

10 s 

1.0 

114 ± 19 

5.9 ± 1.0 


h 2 

2.59t° 0 ° 1 3 X 10 s 

1.0 ± 0.5 x 

10 5 

1.0 

(5060 ± 70 ) 

(260 ± 3) 


Br7 

3-4tS'J X 10 s 

1.9 ± 0.2 x 

10 s 

1.0 

126 ± 13 

6.4 ± 0.7 

MCG+08-11-002 

stars - disk 

> 8.7±|'J x 10 7 

1.0 ± 0.1 X 

10 s 

1.0 ± 0.1 

85 ± 19 

6.1 ± 1.3 


h 2 

6.9ij'g X 10 7 

1.7 ± 0.1 X 

10 s 

0.7 ± 0.1 

61 ± 9 

4.4 ± 0.6 


Br7 

2-3t 1:1 X 10 7 

2.8 ± 0.2 x 

10 s 

0.9 ± 0.1 

24 ± 9 

1.7 ± 0.6 

NGC 2623 

stars - disk 

> 2.9±“;j. x 10 s 

1.4 ± 0.2 x 

10 s 

1.0 ± 0.1 

127 ± 14 

9.1 ± 1.0 

UGC5101 

stars - disk 

> 6.5±l;f x 10 s 

9.4 ± 0.4 x 

10 s 

1.1 ± 0.1 

80 ± 22 

3.0 ± 0.8 

Mrk273 N 

[Fe II] 

1.0 ± 0.1 x 10 9 

1.1 ± 0.1 X 

10 9 

2.0 

89 ± 11 

3.4 ± 0.4 

NGC 6240 N 

stars - disk 

V 

00 

bo 

1 + 
00 

M +1 

X 

0 

00 

1.2 ± 0.1 X 

10 8 

0.63 ± 0.04 

235 ± 8 

13.7 ± 0.5 

NGC 6240 S 

stars - disk 

> 8.7 ± 0.3 x 10 s 

2.8 ± 0.1 X 

10 s 

1.5 ± 0.1 

212 ± 9 

12.3 ± 0.5 

IRASF17207-0014 E 

[Fe II] 

2.5t°o°3 3 x 1q9 

5.3 ± 0.3 X 

10 7 

0.0 ± 0.1 

346 ± 10 

11.6 ± 0.3 

IRASF17207-0014 W 

[Fe II] 

8.9±t:l X 10 7 

6.1 ± 0.4 x 

10 s 

1.1 ± 0.1 

142 ± 9 

4.8 ± 0.3 


a Some fits were run with 7 fixed to match the light profile; these are thus reported with no error bars. 

b Radii in parentheses are artificially large because the mass of the black hole is comparable to or larger than the fitted nuclear disk. 
Presumably the mass profile of the galaxy changes outside of our field of view, which would reduce the radius of the sphere of influence. 
In the case of IIIZw035 H 2 , non-rotational kinematics unseen in other tracers are likely affecting the fit. 
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Fig. 14.— a) Map of Br 7 velocities of the inner region of MCG+08-11-002. b) Map of model velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 15.— a) Map of H 2 velocities of the inner region of MCG+08-11-002. b) Map of model velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 16.— a) Map of stellar velocities of the inner region of NGC2623. b) Map of model stellar velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 17. Left: Map of the symmetrized v rms = [y 2 + cr 2 ) 0 5 of NGC2623. Center: Model v rm s of galaxy 
with black hole. Right: Residual map of model - data divided by errors in v rms . All three panels are centered 
on the black hole and have been rotated so the major axis is horizontal. 
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Fig. 18.— a) Map of stellar velocities of the inner region of UGC5101. b) Map of model stellar velocities of 
best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 19.- Left: Map of the symmetrized v rms = [y 2 + cr 2 ) 0 5 of UGC5101. Center: Model v r ms of galaxy 
with black hole. Right: Residual map of model - data divided by errors in v rms . All three panels are centered 
on the black hole and have been rotated so the major axis is horizontal. 
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Fig. 20.— a) Map of stellar velocities of the inner region of NGC 6240N. b) Map of model stellar velocities 
of best fit black hole model, c) Residual map of best-fit model divided by velocity errors. In panels a) and 
b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
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Fig. 21. - a) Map of [Fe ii] velocities of the inner region of IRAS F17207-0014. b) Map of model velocities 
of best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) 
and b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
We have masked out regions in order to fit only the eastern black hole. 
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Fig. 22.— a) Map of [Fe n] velocities of the inner region of IRAS F17207-0014. b) Map of model velocities 
of best fit black hole model, c) Residual map of best-fit model divided by errors in velocity. In panels a) 
and b), velocity contours are marked in white. In each panel, a red X marks the position of the black hole. 
We have masked out regions in order to fit only the western black hole. 
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